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We compare the results of constructing binary black hole initial data with three different de- 
compositions of the constraint equations of general relativity. For each decomposition we compute 
the initial data using a superposition of two Kerr-Schild black holes to fix the freely specifiable 
data. We find that these initial-data sets differ significantly, with the ADM energy varying by as 
much as 5% of the total mass. We find that all initial-data sets currently used for evolutions might 
contain unphysical gravitational radiation of the order of several percent of the total mass. This 
is comparable to the amount of gravitational-wave energy observed during the evolved collision. 
More astrophysically realistic initial data will require more careful choices of the freely specifiable 
data and boundary conditions for both the metric and extrinsic curvature. However, we find that 
the choice of extrinsic curvature affects the resulting data sets more strongly than the choice of 
conformal metric. 

PACS numbers: 04.25.Dm, 04.70.Bw 



I. INTRODUCTION 

Numerical evolutions of black holes have been im- 
proved slowly but steadily over the last few years and now 
first attempts are being made to extract physical infor- 
mation from these evolutions. Most notably one wants to 
predict the gravitational radiation emitted during black 
hole coalescence [Q, ^, ^ . 

The quality of the initial data will be crucial to the 
success of the predictions of the gravitational wave forms. 
Unphysical gravitational radiation present in the initial 
data will contribute to the gravitational waves computed 
in an evolution and might overwhelm the true gravita- 
tional wave signature of the physical process under con- 
sideration. Therefore an important question is how to 
control the gravitational wave content of initial-data sets, 
and how to specify astrophysically relevant initial data 
with the appropriate gravitational wave content, for e.g. 
two black holes orbiting each other. Unfortunately, as- 
sessing and controlling the gravitational wave content of 
initial-data sets is not well understood at all. 

The mere construction of an initial-data set alone is 
fairly involved, since every initial-data set must satisfy a 
rather complicated set of four partial differential equa- 
tions, the so-called constraint equations of general rel- 
ativity. The question of how to solve these equations, 
and how to specify initial data representing binary black 
holes in particular, has received considerable attention. 
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We consider in this paper three different approaches 
that transform the constraint equations into elliptic 
equations: The conformal transverse-traceless (TT) 
decomposition^ , the physical TT decomposition |^ 
and the conformal thin sandwich decomposition^ . These 
decompositions split the variables on the initial-data sur- 
face into various pieces in such a way that the constraint 
equations determine some of the pieces while not restrict- 
ing the others. After these freely specifiable pieces are 
chosen, the constraint equations are solved and the re- 
sults are combined with the freely specifiable pieces to 
yield a valid initial-data set. 

Any reasonable choice for the freely specifiable pieces 
will lead to a valid initial-data set. Furthermore, any one 
of these decompositions can generate any desired initial- 
data set, given the correct choices of the freely specifiable 
pieces. However, it is not clear what choices of freely 
specifiable pieces lead to initial-data sets with the desired 
properties. 

The decompositions we consider here lead to four 
coupled nonlinear elliptic partial differential equations. 
Since such equations are difficult to solve, the early ap- 
proach to constructing initial data was pragmatic: One 
used the conformal TT decomposition with additional 
restrictions on the freely specifiable pieces, most notably 
conformal flatness and maximal slicing. These assump- 
tions decouple the constraints and allow for analytical 
solutions to the momentum constraints, the so-called 
Bowen-York extrinsic curvature^, [l^. All that re- 
mains is to solve a single elliptic equation, the Hamilto- 
nian constraint. This approach has been used in several 
variations 

However, these numerical simplifications come at a 
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cost. The freely specifiable pieces have been restricted 
to a small subset of all possible choices. One therefore 
can generate only a subset of all possible initial-data sets, 
one that might not contain the desired astrophysically 
relevant initial-data sets. 

Over the last few years there have been additional de- 
velopments: Post-Newtonian results have indicated that 
binary black hole metrics are not conformally flatjl^, p^ . 
With certain restrictions on the slicing, it has also been 
shown that a single stationary spinning black hole can- 
not be represented with a conformally flat spatial met- 
~ 
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In [19 



it was shown that conformally flat 
initial data sets for spinning binary black holes contain 
an unphysical contamination. Moreover, computations 
in spherical symmetry pc|] indicated that initial-data sets 
depend strongly on the choice of the extrinsic curvature 
and that the use of the Bowen-York extrinsic curvature 
might be problematic. 

Therefore it is necessary to move beyond conformally 
flat initial data and to explore different choices for the 
extrinsic curvature. Matzner et al||2l[| proposed a non- 
fiat conformal metric based on the superposition of two 
Kerr-Schild metrics; a solution based on this proposal 
was obtained in [p2| . This work demonstrated the ex- 
istence of solutions to the 3D set of equations, but did 
not examine the data sets in any detail. Refs. ||2^, |2j] 
obtained solutions to a similar set of equations during 
the computation of quasi-circular orbits of binary black 
holes. However, these works assumed conformal flatness. 

In this paper we present a code capable of solving 
the three above-mentioned decompositions of the con- 
straint equations for arbitrary choices of the freely speci- 
fiable pieces. This code is based on spectral methods 
which have been used successfully for several astrophys- 
ical problems (see e.g. @ |2|, |6[|^, |2|, JsO] ) . Our 
code is described in detail in a separate paper 1 31 1. 

We compute solutions of the different decompositions 
for the non-flat conformal metric proposed in Ref. 
Each decomposition has certain choices for the freely 
specifiable pieces and boundary conditions that seem 
"natural" and which we use in our solutions. We com- 
pare the computed initial-data sets with each other and 
with the "standard" conformally-flat solution using the 
Bowcn-York extrinsic curvature. Our major results con- 
firm that 

1. the different decompositions generate different 
physical initial-data sets for seemingly similar 
choices for the freely specifiable pieces. 

2. the choice of extrinsic curvature is critical. 

The first result is certainly not unexpected, but each of 
these factors can cause relative differences of several per 
cent in gauge-invariant quantities like the ADM-energy. 

We also find that the conformal TT/physical TT 
decompositions generate initial-data sets with ADM- 
energies 2 — 3% higher than data sets of the conformal 
thin sandwich decomposition. We demonstrate that this 



higher ADM-energy is related to the choice of the freely 
specifiable part of the extrinsic curvature. In addition, 
we find that the solutions depend significantly on the 
boundary conditions used. 

The paper is organized as follows. In the next sec- 
tion we describe the three decompositions. Section III 
explains how we choose the freely specifiable data within 



IV we describe and test 



each decomposition. In section 
our elliptic solver. Section ^ presents our results, which 
we discuss in section VI. 



II. DECOMPOSITIONS OF EINSTEIN'S 
EQUATIONS AND THE CONSTRAINT 
EQUATIONS 

A. 3+1 Decomposition 

In this paper we use the standard 3-1-1 decomposition 
of Einstein's equations. We foliate the spacetime with 
t — const hypersurfaces and write the four-dimensional 
metric as 

^'^Us^ = -N^ dt^ + 7y {dx' + N'dt){dx^ + N^dt), (1) 

where represents the induced 3-metric on the hyper- 
surfaces, and N and iV' represent the lapse function and 
the shift vector, respectively. We define the extrinsic cur- 
vature Kij on the slice by 



(2) 



where ^^^g is the space-time metric, n the unit normal 
to the hypersurface, and ± denotes the projection op- 
erator into the t = const slice. Einstein's equations di- 
vide into constraint equations, which constrain the data 
(7y , K^^) on each hypersurface, and into evolution equa- 
tions, which determine how the data {jij,K'^^) evolve 
from one hypersurface to the next. The constraint equa- 
tions are 



R + K^ - KijK'^ ^ WnGp 
Vj {K'^ - f^K) = 8nGf. 



(3) 
(4) 



Eq. (^) is called the Hamiltonian constraint, and Eq. (^) 
is referred to as the momentum constraint. K = jijK^^ is 
the trace of the extrinsic curvature, V and R denote the 
three dimensional covariant derivative operator and the 
Ricci scalar compatible with 7^. p and are the energy 
and momentum density, respectively. Both vanish for the 
vacuum spacetimes considered here. 
The evolution equation for 7^ is 



(5) 



which follows from Eq. (|^). There is a similar albeit 
longer equation for dfKij which we will not need in this 
paper. The choices of N and N'^ are arbitrary. One can 
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in principle use any lapse and shift in the evolution off 
the initial-data surface, although some choices of lapse 
and shift are better suited to numerical implementation 
than others. 

Later in this paper we will often refer to the trace-free 
piece of Eq. (H). Denote the tracefree piece of a tensor 
by TF(.), and define 7 = det7y . From Eq. (^) and the 
fact that 6lnj = ^'^''S-fki, it follows that 

TF(da,j) = -f^/^dt (7~'/'7y) = + (LTV),,. 

(6) 

Here Aij — K^j — ■^JijK denotes the trace- free extrinsic 
curvature, and 

{LNy^ = S/'N^ + V^N' - ^7*^ VfciV^ (7) 

L always acts on a vector, so the 'N' in (LA^)*-' denotes 
the shift vector A^* and not the lapse TV. 



B. Decomposition of tlie constraint equations 



1. Conformal TT Decomposition 

In this decomposition one first conformally transforms 
the traceless extrinsic curvature, 

yiy=^-ioiy\ (10) 

and then applies a TT decomposition with respect to the 
background metric ^ij : 

A'^^ = A^^rj. + {tXyi . (11) 

The operator L is defined by Eq. (|^) but using the con- 
formal metric 7^, and derivatives associated with 7^. 
Ai^rp is transverse with respect to the conformal metric, 
VjAjlrp = 0, and is traceless. 

Substituting Eqs. (^) and ( [III ) into the momentum 
constraint (^), one finds that it reduces to an elliptic 
equation for X^, whereas Al^j, is unconstrained. 

In order to specify the transverse-traceless tensor A'^j, 
one usually has to construct it from a general symmet- 
ric trace-free tensor M*-' by subtracting the longitudi- 
nal piece. As described in |32|] one can incorporate the 
construction of Al^rp from M*-' into the momentum con- 
straint, arriving at the following equations: 



Equations (||) and (^ constrain four degrees of free- 
dom of the 12 quantities (7^,/^'^). However, it is not 
immediately clear which pieces of 7^ and K''^ are con- 
strained and which pieces can be chosen at will. Several 
decompositions have been developed to divide the 12 de- 
grees of freedom into freely specifiable and constrained 
pieces. We will now review some properties of the three 
decompositions we consider in this paper. 

All three decompositions follow the York-Lichnerowicz 
approach and use a conformal transformation on the 
physical 3- metric jij , 

jij =ip%j. (8) 

-0 is called the conformal factor, jij the background met- 
ric or conformal metric. We will denote all conformal 
quantities with a tilde. In particular, V is the covariant 
derivative operator associated with 7^ , and Rij and R 
are the Ricci tensor and Ricci scalar of 7^. 

The extrinsic curvature is split into its trace and trace- 
free part, 

K'^ = A'^ + -j^^K. (9) 
3 

The three decompositions of the constraint equations we 
discuss in this paper differ in how A'^^ is decomposed. For 
each decomposition, we discuss next the relevant equa- 
tions, and describe how we choose the quantities one has 
to specify before solving the equations. We use the con- 
ventions of p^ . 



o 1^ o 

(12) 

- ^V^V'if + VjlvP' = 8TTGilj^°j\ (13) 

where A*-' and the operator Al are defined by 

A'^ = {Lvy^ + M'^ (14) 

and 

Ai^V = ^,{lvy\ (15) 

After solving these equations for ip and V^, one ob- 
tains the physical metric jij from (^ and the extrinsic 
curvature from 

K'^ = V"^"^'^' + -ip'^f^K. (16) 
3 

We will refer to Eqs. (|l|) and (|l|) together with (|l|), 
( |l6| ) and (||) as the conformal TT equations. In these 
equations we are free to specify the background metric 
7ij, the trace of the extrinsic curvature K, and a sym- 
metric traceless tensor JVP^ . The solution will contain 
a contribution that removes the longitudinal piece from 
M*-' and the piece that solves the momentum constraint 
if M*-' were transverse-traceless. 

This decomposition has been the most important in 
the past, since if one chooses a constant K and if one 
considers vacuum spacetimes then the momentum con- 
straint (nM decouples from the Hamiltonian constraint 
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([14) ■ Moreover, if one assumes conformal flatness and 
AP^ = 0, it is possible to write down analytic solutions 
to Eq. ([l3|), the so-called Bowen-York extrinsic curva- 
ture. In that case one has to deal with only one elliptic 
equation for tjj. The Bowen-York extrinsic curvature can 
represent multiple black holes with arbitrary momenta 
and spins. One can fix boundary conditions for ^ by re- 
quiring that the initial-data slice be inversion symmetric 
at both throats In that case one has to modify 

the extrinsic curvature using a method of images. We will 
include initial-data sets obtained with this approach be- 
low, where we will refer to them as inversion symmetric 
initial data. 

Reasonable choices for the freely specifiable pieces 7^ , 
K, M"^^ will lead to an initial-data set {'-fij,K'^^) that 
satisfies the constraint equations. How should we choose 
all these functions in order to obtain a desired physical 
configuration, say a binary black hole with given linear 
momenta and spins for the individual holes? We can gain 
insight into this question by considering how the confor- 
mal TT decompositions can recover a known solution. 

Suppose we have a known solution (70 ij , Kq ) of the 
constraint equations. Denote the trace and trace-free 
parts of this extrinsic curvature by Kq and Aq"* , respec- 
tively. If we set 



lij = 70 i] , 



^0 



then 



= 1, = 



(17) 



(18) 



trivially solve Eqs. (|12- 13). Note that we have to set M*-' 
equal to the trace-free part of the extrinsic curvature. 

Now suppose we have a guess for a metric and an ex- 
trinsic curvature, which — most likely — will not satisfy 
the constraint equations @} and . Set 7ij to the guess 
for the metric, and set K and M^^ to the trace and 
trace-free piece of the guess of the extrinsic curvature. 
By solving the conformal TT equations we can compute 
(7y , X*-' ) that satisfy the constraint equations. If our 
initial guess is "close" to a true solution, we will have 
■0 « 1 and « 0, so that 7^ and K^^ will be close to 
our initial guess. 

Thus one can guess a metric and extrinsic curvature as 
well as possible and then solve the conformal TT equa- 
tions to obtain corrected quantities that satisfy the con- 
straint equations. 

An artifact of the conformal TT decomposition is that 
one has no direct handle on the transverse traceless piece 
with respect to the physical metric. For any vector X*, 



Thus, Eqs. (|l|) and (|T|) imply 



For any symmetric traceless tensor S^^ 



(19) 



(20) 



(21) 



Therefore the first term on the right hand side of Eq. ( |2(]| ) 
is transverse-traceless with respect to the physical metric, 



7, 



(22) 



However, the second term on the right hand side of 
Eq. (|2^) is conformally weighted. Therefore, Eq. (|20| ) 
does not represent the usual TT decomposition. 



2. Physical TT Decomposition 

In this case one decomposes the physical traceless ex- 
trinsic curvature directly: 



(23) 



As above in the conformal TT decomposition, the mo- 
mentum constraint becomes an elliptic equation for AT*. 
We can again incorporate the construction of the sym- 
metric transverse traceless tensor from a general 
symmetric tensor M*^ into the momentum constraint. 
Then one obtains the physical TT equations: 



(24) 



2~ 



4-i 



(25) 



where A''^ is defined by 

A'^ = {Lvy^ +^p'^M'^. 



(26) 



When we have solved (24) and ( p5[ ) for ^ and V'', the 
physical metric is given by (^, and the extrinsic curva- 
ture is 



A'- 



(27) 



We are free to specify the background metric 7^, the 
trace of the extrinsic curvature and a symmetric trace- 
less tensor M*-' . As with the conformal TT equations, the 
solution will contain a contribution that removes the 
longitudinal piece from AP^ and a piece that solves the 
momentum constraint if M^^ were transverse-traceless. 

These equations can be used in the same way as the 
conformal TT equations. Guess a metric and extrinsic 
curvature, set 7^ to the guess for the metric, and K 
and M'^ to the trace and trace-free pieces of the guess 
for the extrinsic curvature. Then solve the physical TT 
equations to obtain a corrected metric 7^ and a corrected 
extrinsic curvature K^^ that satisfy the constraint equa- 
tions. 

The transverse traceless piece of K^^ (with respect to 
7ij) will be the transverse traceless piece of ■(/'"^'^M'^. 
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One can also easily rewrite the physical TT equations 
such that il}~^^M'^^ can be freely chosen instead of M^^ . 
So, in this decomposition we can directly control the TT 
piece of the physical extrinsic curvature. We have chosen 
to follow since it seems somewhat more natural to 
specify two conformal quantities, 7^ and M*-' than to 
specify one conformal and one physical quantity. 



3. Conformal thin sandwich decomposition 

The conformal and physical TT decompositions rely 
on a tensor splitting to decompose the trace-free part of 
the extrinsic curvature. In contrast, the conformal thin 
sandwich decomposition simply defines A^^ by Eq. ( |lO| ) 
and the decomposition 



(28) 



where w*-' is symmetric and tracefree. Eq. (^8|) is moti- 
vated by Eq. (|): If one evolves an initial-data set with 
A^^ of the form (pSl) using as lapse and shift 



N 



then 



(29) 



(30) 



Therefore, the decomposition ( pSD is closely related to the 
kinematical quantities in an evolution. Although d and 
/?' are introduced in the context of initial data, one usu- 
ally refers to them as the "conformal lapse" and "shift" . 
While the form of Eq. ( psf ) is similar in form to the con- 
formal and physical TT decompositions, there are differ- 
ences. In particular, ■0*-' is not divergenceless. 

Within the conformal thin sandwich decomposition, 
the constraint equations take the form: 



(31) 



aVj -u'^ = 16TTG6i'i(j'"f (32) 



Having solved Eqs. (^) and (|3^) for ip and the vector /3% 
one obtains the physical metric from (^) and the extrinsic 
curvature from 
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(33) 



In this decomposition we are free to specify a conformal 
metric 7^, the trace of the extrinsic curvature K, a sym- 
metric trace-free tensor u^^ and a function d. 

It seems that the conformal thin sandwich decomposi- 
tion contains additional degrees of freedom in the form of 



the function d and three additional unconstrained com- 
ponents of u^^ . This is not the case. The longitudinal 
piece of u^^ corresponds to the gauge choice of the ac- 



tual shift vector used in an evolution. Thus 



really 



only contributes two degrees of freedom, just like M*-' 
in the conformal and physical TT decompositions. Fur- 
thermore, we can reach any reasonable physical solution 
(7y , K^^) with any reasonable choice of a; each choice of 
d simply defines a new decomposition. A forthcoming ar- 
ticle by YorkjS^ will elaborate on these ideas. Note that 
for a = 1/2 we recover the conformal TT decomposition. 

Let us now turn to the question of how one should 
pick the freely specifiable data in the conformal thin 
sandwich approach. We motivate our prescription again 
by considering how to recover a known spacetime: As- 
sume we are given a full four-dimensional spacetime with 
3-1-1 quantities 70 ij, K^^ , iVg and A^o- Further assume 
the spacetime is stationary and the slicing is such that 
dt"fij — dtKij — 0. An example for such a situation is a 
Kerr black hole in Kerr-Schild or Boyer-Lindquist coor- 
dinates. 



Using i9t7oy = in Eq. ^ 
trace-free extrinsic curvature 



yields a relation for the 



(34) 



This is a decomposition of the form ( p8[ ) with u^^ = 0. 
Therefore, if we choose the freely specifiable data for the 
conformal thin sandwich equations as 



K = Ko, ii'^ = 0, 



(35) 



and if we use appropriate boundary conditions, then the 
solution of the conformal thin sandwich equations will 
be -0 = 1 and — Nq. As part of the solution, we ob- 
tain the shift vector needed for an evolution to produce 
TF(9t7ij) = 0. Not needing a guess for the trace-free 
extrinsic curvature, and having the solution automat- 
ically provide an initial shift for evolution, make the con- 
formal thin sandwich equations very attractive. 

In order to generate initial-data slices that permit an 
evolution with zero time derivative of the conformal met- 
ric — a highly desirable feature for quasi-equilibrium 
data, or for a situation with holes momentarily at rest — 
one can proceed as follows: Set %j and K to the guess for 
the metric and trace of extrinsic curvature, respectively. 
Set d to the lapse function that one is going to use in 
the evolution, and set u^^ — 0. If these guesses are good, 
the conformal factor -ip will be close to 1, and N — ip^d 
as well as A^' — fi^ give us the actual lapse function and 
shift vector to use in the evolution. 
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III. CHOICES FOR THE FREELY 
SPECIFIABLE DATA 

A. Kerr-Schild coordinates 

We base our choice for the freely specifiable data on a 
superposition of two Kerr black holes in Kerr-Schild co- 
ordinates. In this section we describe the Kerr-Schild so- 
lution and collect necessary equations. We also describe 
how we compute the 3-metric, lapse, shift and extrinsic 
curvature for a boosted black hole with arbitrary spin. 

A Kerr-Schild metric is given by 



(36) 



where 77^^ is the Minkowski metric, and is a null- vector 
with respect to both the full metric and the Minkowski 
metric: g^'^l^lp = rj^^'^l^jli, = 0. The 3-metric, lapse and 
shift are 



velocity v'^. However, if the evolution uses the shift vector 
TV* -f f% the coordinates will move with the black hole, 
and the hole will be at rest in coordinate space. The 
spacetime is nonetheless different from a Kerr black hole 
at rest. The ADM-momentum will be P 



ADM 



where M is the rest- mass of the hole and 7 = (1- 



-V 



-fMv\ 
•2^-1/2^ 



B. Freely specifiable pieces 

We want to generate initial data for a spacetime con- 
taining two black holes with masses Ma.b^ velocities va,b 
and spins Mao, a and Mbolb- 

We follow the proposal of Matzner et al j2^, and 
base our choices for the freely specifiable choices on two 
Kerr-Schild coordinate systems describing two individual 
black holes. The first black hole with label A has an 
associated Kerr-Schild coordinate system with metric 



7u 
N ■ 



2Hlilj, 



TV* = - 



l + 2HlHt' 



(37) 
(38) 

(39) 



For a black hole at rest at the origin with mass M and 
angular momentum Ma, one has 



H = 



irest 



'rest 



+ (a • x)^ ' 

(1, ^rest), 

rx — a X X + {a ■ x)a/r 



j,2 _j_ g2 



with 



+ (a • f)' 



1/2 



(40) 
(41) 
(42) 

(43) 



For a nonrotating black hole with a — 0, H has a pole 
at the origin, whereas for rotating black holes, r has a 
ring singularity. We will therefore have to excise from 
the computational domain a region close to the center of 
the Kerr-Schild black hole. 

Under a boost, a Kerr-Schild coordinate system trans- 
forms into a Kerr-Schild coordinate system. Applying 
a Lorentz transformation with boost velocity to ^jj'^^*, 
we obtain the null-vector Lof the boosted Kerr-Schild 
coordinate system. Eqs. ( p7[ - |39| ) give then the boosted 
3-metric, lapse, and shift. Since all time-dependence is 
in the uniform motion, evolution with lapse TV and shift 
N'- yields dt"fij = —v^dk^ij, and from Eq. (||) one can 
compute the extrinsic curvature 



1 

2N 



{v''^kl^J + V.TV, -I- V.TV,) 



(44) 



If this initial-data set is evolved with the shift TV*, the 
black hole will move through the coordinate space with 



lAij = Sij + 2Ha Iai I 



(45) 



and with an extrinsic curvature KAij, a lapse Na and a 
shift N\. The trace of the extrinsic curvature is Ka- All 
these quantiti es can be computed as described in the pre- 
vious section. III A. The second black hole has a similar 



set of associated quantities which are labeled with the 
letter B. 

For all three decompositions, we need to choose a con- 
formal metric and the trace of the extrinsic curvature. 
We choose 

= + 2Ha I a ^lA3+ 2Hb Ib.Ib j (46) 
K^Ka + Kb (47) 

The metric is singular at the center of each hole. There- 
fore we have to excise spheres around the center of each 
hole from the computational domain. We now specify 
for each decomposition the remaining freely specifiable 
pieces and boundary conditions. 



1. Conformal TT and physical TT decompositions 

For the conformal TT and physical TT decompositions 
we will be solving for a correction to our guesses. As guess 
for the trace-free extrinsic curvature, we use a superpo- 
sition 



K 



Ak ' 



K 



Bk 



1 



6'i{KA + Kb, 



(48) 



TVf*-' is symmetric and trace-free with respect to the con- 
formal metric, "fijM^^ = 0. Solving for a correction only, 
we expect that ip ^ 1 and 0, hence we use Dirichlet 
boundary conditions 



^ = 1, 



V = 0. 



(49) 
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2. Confomnal thin sandwich 

For conformal thin sandwich, we restrict the discussion 
to either two black holes at rest, or in a quasi-circular or- 
bit in corotating coordinates. In these cases, one expects 
small or even vanishing time-derivatives, dt^iij ~ 0, and 
so Eq. (pG) yields the simple choice 



0. 



(50) 



The conformal 3-mctric and the trace of the extrinsic 
curvature are still given by Eqs. ( ^ ) and (^). Orbit- 
ing black holes in a corotating frame will not move in 
coordinate space, therefore we do not boost the individ- 



ual Kerr-Schild metrics in this decomposition: v 



A/B 



0. 



The lapse functions Na/b and the shifts N^^g are also 



for unboosted Kerr-Schild black holes. 
We use Dirichlet boundary conditions: 





= 1 


all boundaries 


(51a) 




-NX 


sphere inside hole A 


(51b) 




^Nh 


sphere inside hole B 


(51c) 




= fl X r 


outer boundary 


(51d) 



Eq. (51d) ensures that we are in a corotating reference 
frame; the cross-product is performed in flat space, and 
r2 = corresponds to two black holes at rest. Close to 
the holes we force the shift to be the shift of a single black 
hole in the hope that this choice will produce a hole that 
is at rest in coordinate space. 

For the conformal lapse we use 



a = NA + NE 



1 



or 



a 



Na Nb. 



(52) 



(53) 



The first choice of a follows the philosophy of adding 
quantities of each individual hole. However, a of Eq. ( |5^ ) 
becomes negative sufficiently close to the center of each 
hole and is therefore a bad choice if the excised spheres 
are small. The choice (53) does not change sign and has 
at large distances the same behavior (same 1 /r term) as 



IV. NUMERICAL IMPLEMENTATION 



We implemented an 
three decompositions 
generality. The solver 
can handle nontrivial 
dospectral collocation, 
in each subdomain as 
This elliptic solver is 
paper 



elliptic solver that can solve all 
we described above in complete 
uses domain decomposition and 
topologies. It is based on pseu- 
that is, it expresses the solution 
an expansion in basis functions, 
described in detail in a separate 




FIG. 1: Structure of domains. Spherical shells around each 
excised sphere are surrounded by 43 rectangular blocks and 
another spherical shell. The rectangular blocks touch each 
other and overlap with all three spherical shells. 



From the computational domain we excise two spheres 
containing the singularities of the Kerr-Schild metric 
close to the center of each hole. Around each of the 
excised spheres we place a spherical shell. These shells 
are patched together with 5 x 3 x 3 = 45 rectangular 
blocks, with the two blocks at the location of the spheres 
removed. Around these 43 blocks, another spherical shell 
is placed that extends far out, typically to an outer radius 
of lO^Af. In the rectangular blocks, we expand in Cheby- 
shev polynomials, while in the spheres we use Chebyshev 
polynomials radially and spherical harmonics for the an- 
gular variables. This setup is depicted in Fig. |l|. 

The domain decomposition in Fig. ^ is fairly compli- 
cated. Even if the shells were made as large as possible, 
they do not cover the full computational domain when 
the excised spheres are close together. Thus additional 
subdomains are needed in any case. Choosing the 43 
cubes as depicted allows for relatively small inner shells 
and for a relatively large inner radius of the outer shell. 
Thus each shell covers a region of the computational do- 
main in which the angular variations of the solution are 
fairly low, allowing for comparatively few angular basis- 
functions. 

The code can handle a general conformal metric. In 
principle, the user needs to specify only 7.^ . Then the 
code computes 7'^ , and — using numerical derivatives — 
the Christoffel symbols, Ricci tensor and Riemann scalar. 
For the special case of the Kerr-Schild metric of a single 
black hole and the superposition of two Kerr-Schild met- 
rics, Eq. (^), we compute first derivatives analytically 
and use numerical derivatives only to compute the Rie- 
mann tensor. 

The solver implements Eqs. (|l^) and (^3|) for the con- 
formal TT decomposition, Eqs. ( p^ ) and ( |25| ) for the 
physical TT decomposition, and Eqs. (31) and (p^) for 
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the conformal thin sandwich decomposition. 

After solving for (■0, V^) [conformal TT and physical 
TT], or {ip,f3^) [thin sandwich] we compute the physical 
metric 7ij and the physical extrinsic curvature K^^ of the 
solution. Utilizing these physical quantities (7^, K^^), we 
implement several analysis tools. We evaluate the con- 
straints in the form of Eq. (|^) and we compute ADM- 
quantities and we search for apparent horizons. Note 
that these analysis tools are completely independent of 
the particular decomposition; they rely only on 7^ and 

Next we present tests ensuring that the various systems 
of equations are solved correctly. We also include tests of 
the analysis tools showing that we can indeed compute 
constraints, ADM-quantities and apparent horizons with 
good accuracy. 



A. Testing the conformal TT and physical TT 
decompositions 

We can test the solver by conformally distorting a 
known solution of the constraint equations. Given a solu- 
tion to the constraint equations (70 ij , K^'' ) pick functions 



* > 0, 



and set 



and 



7y 
K ^Ko, 



(54) 



(55) 
(56) 



^0 



for conformal TT or 



jI'Ko] ~^\hoWy^ (57) 



(58) 



for physical TT. With these freely specifiable pieces and 
appropriate boundary conditions, a solution of the con- 
formal TT equations (|l|), (|l|) or the physical TT equa- 
tions (|23), (M) will be 



V" = * 



(59) 
(60) 



From Eq. (Q) we recover our initial metric 70 ij , and from 
Eq. dH) [conformal TT] or Eq. (|^) [physical TT] we 



recover the extrinsic curvature Kq . 

In our tests we used the particular choices 

8(r ^ 2) 



^ = 1 



36 + a;2-H0.V + l.3(z-l) 



W = 



50(7 



(64 + r4) 



i-y,x, 1). 



(61) 
(62) 




FIG. 2: Plot of the functions * and from Eqs. ^) and 
(p2[). The solid line depicts ^I* along the positive x-axis, the 
short dashed line depicts along the negative z-axis. The 
long dashed line is a plot of along the positive x-axis. 



These functions are plotted in Fig. ||. 5* varies between 
0.8 and 1.5, varies between ±0.5, and both take their 
maximum values around distance 7 from the center 
of the hole. We used for (70 ij , -ft^o'' ) a single, boosted, 
spinning black hole in Kerr-Schild coordinates. 

Figure ^ shows results of testing the conformal TT 
decomposition on a single spherical shell. The numer- 
ical solution {ip, V) converges to the analytic solutions 
(^', W'') exponentially with the number of basis functions 
as expected for a properly constructed spectral method. 
Moreover, the reconstructed metric and extrinsic curva- 
ture satisfy the constraints. 

Now we test the solver for the physical TT decompo- 
sition, and demonstrate that we can correctly deal with 
multiple domains. In this example the computational do- 
main is covered by an inner spherical shell extending for 
1.5M < r < lOM. This shell is surrounded by 26 rectan- 
gular blocks that overlap with the shell and extend out 
to x,y,z = ±25M. Finally another spherical shell covers 
the region 20Af < r < lO'^M. As can be seen in Fig. ||, 
the solution converges again exponentially. 

For realistic cases we do not know the analytic solution 
and therefore need a measure of the error. Our major tool 
will be the change in results between different resolution. 
In particular we consider the L2 norm of the point-wise 
differences of the solution at some resolution and at the 
next lower resolution. This diagnostic is labeled by cir- 
cles in Fig. |4[ Since the solution converges exponentially, 
these circles essentially give the error of the lower of the 
two resolutions. 

In addition to testing the equations, this example tests 
domain decomposition and the integration routines for 
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FIG. 3: Testin g th e solver for the conformal TT decom- 
position. Eqs. (|l^p^ with freely specifiable data given 
by Eqs. (pq-M) are solved in a single spherical shell with 
1.5M < r < lOM. N is the cube root of the total number 
of unknowns. Plotted are the L2-norms oi 4> — 'i' , V — W, 
and the residuals of Hamiltonian and momentum constraints, 
Eqs. (|) and (|). 
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FIG. 4: Physical TT decomposition with domain decom- 
position. Eqs. ([24[ ). ( p5| ) with freely specifiable data given 
by Eq. ^ [s^ ) are solved in multiple domains (one in- 
ner spherical shell, 26 rectangular blocks, one outer spherical 
shell). A'' is the cube root of the total number of grid-points. 
dijj denotes the L2-norm of the difference of the solution and 
the solution at next lower resolution. Triangles denote the 
L2-norm of the difference to the analytic solution. The re- 
maining symbols denote the errors of numerically extracted 
AD M- quant it ies . 




the ADM quantities. The ADM quantities are computed 
by the standard integrals at infinity in Cartesian coordi- 
nates. 



E 



ADM 



1 

16^ 



hijj - l3J,i) d^^i^ (63) 

•^(0-^ / {K^' -f'K)S,,SS... (64) 
For the x-component of the linear ADM-momentum, 



^ = e^; in Eq. (|64|). The choice ^ = xcy — ye^ yields 
the z-component of the ADM-like angular momentum as 
defined by York Q . Since the space is asymptotically flat 
there is no distinction between upper and lower indices 
in Eqs. ^ and (|^. Note that Eq. ^ reduces to the 
familiar monopole term 



1 

'2^ 



(65) 



only for quasi-isotropic coordinates. Our outer domain is 
large, but since it does not extend to infinity, we extrap- 
olate r — > CX3. 

For a Kerr black hole with mass M and spin a, that is 
boosted to velocity u, the ADM-quantities will be 



Eadm = iM, 
Padm = ^Mv, 

Jadm — 



11 ^ 



M, 



(66) 
(67) 

(68) 



where 7 = (1 — i;^)~^" denotes the Lorentz factor. 
Eq. ( |68| ) reflects the fact that under a boost, the com- 
ponent of the angular momentum perpendicular to the 
boost-direction is multiplied by 7. 

The example in Fig. || uses v = (0.2,0.3,0.4), and 
a = (-1/4, 1/4, 1/6)M. The evaluation of the angular 
momentum Jz seems to be less accurate since our cur- 
rent procedure to extrapolate to infinity magnifies round- 
off. We plan to improve this in a future version of the 
code. Until then we seem to be limited to an accuracy of 
- 10-6. 



B. Testing conformal thin sandwich equations 

The previous decompositions could be tested with a 
conformally distorted known solution. In order to test 
the conformal thin sandwich decomposition we need to 
find an analytic decomposition of the form (p8|). To do 
this, we start with a stationary solution of Einstein's 
equations and boost it with uniform velocity u*. Denote 
the metric, extrinsic curvature, lapse and shift of this 
boosted spacetime by 70 ij, = + Kq, Nq and 
iV( 

will not find dtjij 
However, all time-dependence of this spacetime is due to 
the uniform motion, so in the comoving reference frame 



0, respectively. Since we boosted the static solution, we 
if we evolve it with the shift Nq. 
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specified by the shift TVq + 
this case, Eq. (||) yields 



we will find dt'^ij = 0. In 



(69) 



If we choose a — Nq and u'^ = 0, the thin sandwich 
equations ( |3l| ) and ( ^ ) will thus be solved by V' = 1 and 
/3* = TVq + u* . Similar to the conformal TT and physical 
TT decomposition above, we can also conformally distort 
the metric 70 ij- Furthermore, we can consider nonvan- 
ishing ■u*-' . We arrive at the following method to test the 
solver for the conformal thin sandwich equations: 

Given a boosted version of a stationary solution with 
shift Nq, lapse Nq, 3-metric 70 ij, trace of extrinsic cur- 
vature Kq, and boost- velocity t;'. Pick any functions 

* > (70) 

W (71) 

and set 

7y■ = ^'-S^J (72) 

K =K^ (73) 
a = ^-^Nq (74) 
u'^ = *4(LoW^)'J (75) 

Then a solution to the thin sandwich equations (|3l[]3^ ) 
will be 

= (76) 
j3' =Nl,+v' + W' (77) 

assuming boundary conditions respecting this solution. 

Figure || shows results of this test for a single spheri- 
cal shell and a Kerr black hole with v = (0.2, —0.3, 0.1), 
a = (0.4, 0.3, 0.1)M. The solution converges to the ex- 
pected analytical result exponentially. In addition, ap- 
parent horizon searches were performed. For the numeri- 
cally found apparent horizons, the apparent horizon area 
Aah as well as the apparent horizon mass 



Mah = 



A 



AH 



167r 



(78) 



were computed. The figure compares Mah to the ex- 
pected value 



(79) 



As described in the apparent horizon finder ex- 

pands the apparent horizon surface in spherical harmon- 
ics up to a fixed order L. For fixed L, the error in the 
apparent horizon mass is dominated by discretization er- 
ror of the elliptic solver at low resolution N . As N is 
increased, the discretization error of the elliptic solver 
falls below the error due to finite L. Then the error in 
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FIG. 5: Testing thin sandwich decomposition witli apparent 
horizon searches. Equations (Blh and (B2|) with freely specifi- 
able data given by Eqs. (|7^)-([75|) are solved in a single spher- 
ical shell. TV and diff as in Fig~M. Apparent horizon searches 
with different surface expansion order L were performed, and 
the errors of the apparent horizon mass Mah are plotted. 



Mah becomes independent of A'^. Since the expansion in 
spherical harmonics is spectral, the achievable resolution 
increases exponentially with L. Note that for exponential 
convergence it is necessary to position the rays in the ap- 
parent horizon finder at the abscissas of Gauss-Legendre 
integration. 



C. Convergence of binary black hole solutions 

Figure ^ present the convergence of the solver in the 
binary black hole case. In this particular example, the 
conformal TT equations were solved for two black holes 
at rest with coordinate separation of lOM. The com- 
putational domain is structured as in Fig. g. The ex- 
cised spheres have radius r^xc = 2M, the inner spherical 
shells extend to radius AM . The rectangular blocks cover 
space up to x,y, z = ±25M, and the outer spherical shell 
extending from inner radius 20M to an outer radius of 
R = lO^M. 

We do not use fall-off boundary conditions at the outer 
boundary; we simply set ip = 1 there. This limits the 
computations presented in this paper to an accuracy of 
order 1/i? ~ 10""^. Figure^ shows that even for the next 
to highest resolution (A^ « 80) the solution will be limited 
by the outer boundary condition. All results presented in 
the following section are obtained at resolutions around 
N « 80. If the need arises to obtain solutions with higher 
accuracy, one can easily change to a fall-off boundary 
condition, or just move the outer boundary further out. 
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FIG. 6: Binary black hole with conformal TT decomposition. 
The residuals of several quantities are plotted as a function 
of the cube root of the total number of grid points, diff as 
in Fig. ^, Ham and Mom are the residuals of Hamiltonian 
and momentum constraints. Eadm denotes the difference 
between ADM-energy at resolution A'^ and ADM-energy at 
highest resolution. 



FIG. 7: The conformal factor ij) along the axis connecting 
the holes for several decompositions, x measures the distance 
from the center of mass, so that the excised sphere is located 
between 3 < a; < 7. mConlTT is explained below in section 



V C . The solution of PhysTT is not plotted since it is within 
the line thickness of ConfTT. The insert shows an enlargement 
for large x. 



RESULTS 



The purpose of this paper is to compare the initial-data 
sets generated by different decompositions using simple 
choices for the freely specifiable pieces in each decompo- 
sition. We solve 

• ConfTT: Conformal TT equations (|l|) and (|l|) 
with freely specifiable pieces and boundary conditions 
given by Eqs. 1^, (g|) and ® . 

• PhysTT: Physical TT equations (H) and (||) with 
same freely specifiable pieces and boundary conditions as 
ConfTT. 



• CTS: Conformal thin sandwich equations (31) and 
( p^ ) with freely specifiable pieces and boundary condi- 
tions given by Eqs. (|4^), (F^ , ( pO| ) and (Fll) . The lapse 
a is given by either Eq. (p2|)^ or by Eq. (p3|). 

We will apply the terms "ConfTT", "PhysTT" and 
"CTS" only to these particular choices of decomposition, 
freely specifiable pieces and boundary conditions. When 
referring to different freely specifiable pieces, or a decom- 
position in general, we will not use these shortcuts. If we 
need to distinguish between the two choices of a in CTS, 
we will use "CTS-add" for the additive lapse Eq. (H) 
and "CTS-mu lt" fo r the multiplicative lapse Eq. (^3|) . Be- 
low in section V C we will also introduce as a forth term 
"mConfTT" . 



A. Binary black hole at rest 

We first examine the simplest possible configuration: 
Two black holes at rest with equal mass, zero spin, and 
with some fixed proper separation between the apparent 
horizons of the holes. We solve 

• ConfTT 

• PhysTT 

• CTS (with both choices of a). 

In the comparisons, we also include inversion symmetric 
conformally fiat initial data obtained with the conformal- 
imaging formalism. 

We excised spheres with radius r^xc — 2M , which is 
close to the event horizon for an individual Eddington- 
Finkelstein black hole. This results in the boundary con- 
ditions being imposed close to, but within the apparent 
horizons of the black holes. The centers of the excised 
spheres have a coordinate separation oi d — lOM. 

We now discuss the solutions. The conformal factor 
is very close to 1 for each of the three decompositions. 
It deviates from 1 by less than 0.02, indicating that a 
conformal metric based of a superposition of Kerr-Schild 
metrics does not deviate far from the constraint surface. 

Figure ^ presents a plot of the conformal factor along 
the axis through the centers of the holes. One sees that 
■0 is close to 1; however, between the holes ConfTT and 
CTS force ■0 in opposite directions. For CTS, > 1 
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FIG. 8: Black holes at rest: Contour plots of the conformal 
factor V for ConfTT (left) and CTS-add (right). The circles 
denote the excised spheres of radius 2. 



between the holes, for ConfTT, ijj < 1\ The contour 
plots in Fig. || also show this striking difference between 
the decompositions. 

The result of PhysTT was found to be almost iden- 
tical with ConfTT. This is reasonable, since these two 
decompositions differ only in that in one case the TT de- 
composition is with respect to the conformal metric, and 
in the other case the TT decomposition it is with respect 
to the physical metric. Since ■0 « 1, the conformal met- 
ric is almost identical to the physical metric, and only 
minor differences arise. In the following we will often use 
ConfTT/PhysTT when referring to both data sets. 

We performed apparent horizon searches for these 
cases. For all data sets, the apparent horizon is outside 
the sphere with radius 2M, that is outside the coordinate 
location for the apparent horizon in a single hole space- 
time. For ConfTT/PhysTT, the radius of the apparent 
horizon surface is « 2.05M, for CTS it is « 2.15M. We 
computed the apparent horizon area Aah, the apparent 
horizon mass 



Ma 



H 



A 



AH 



IGtt 



(80) 



of either hole, and the combined mass of both holes, 



TO = 2Mah- 



(81) 



There is no rigorous definition of the mass of an indi- 
vidual black hole in a binary black hole spacetime, and 
Eq. ( ^ ) represents the true mass on an individual black 
hole only in the limit of wide separation of the black holes. 
A hard upper limit on the possible gravitational radia- 
tion emitted to infinity during the coalescence process of 
a binary will be 



Empbc = Ea 



DM 



2Aah 
167r 



(82) 



where 2Aah is the combined apparent horizon area of 
both holes. Thus, Ej^ipj^c represents the maximum pos- 
sible radiation content (MPRC) of the initial data. This, 



of course, makes the unlikely assumption that the binary 
radiates away all of its angular momentum. 

We also compute the proper separation £ between the 
apparent horizon surfaces along the straight line con- 
necting the centers of the excised spheres. In order to 
compare different data sets we consider the dimension- 
less quantities 1/Eadm, EADM/m and Emprc/Eadm- 
We will also use Eb/n which will be defined shortly. 

Table || summarizes these quantities for all three de- 
compositions. It also includes results for inversion sym- 
metric initial data, which for black holes at rest reduces 
to the Misner data|^^. The results in Table | are in- 
tended to represent nearly the same physical configura- 
tion. 

From Table || one finds that the black holes have 
roughly the same dimensionless proper separation. How- 
ever, the scaled ADM-energy EADM/m differs by 
as much as 4.7% between the different data sets. 
Emprc/Eadm, which does not depend on any notion 
of individual black hole masses at all, differs by 16% be- 
tween the different data sets. 

The inversion symmetric data has lowest EADM/fn 
and Emprc/Eadm, CTS has somewhat larger values, 
and ConfTT/PhysTT lead to the biggest values for 
EADM/m and Emprc/Eadm ■ This indicates that, rela- 
tive to the sizes of the black holes, ConfTT/PhysTT and 
CTS probably contain some excess energy. 

A slightly different argument uses the binding energy 
which is defined as 



Eb ^ Eadm - 2Mah 



(83) 



where — Mah/^ is the reduced mass. Two Newtonian 
point masses at rest satisfy 



E^ 



1 



t/r 



(84) 



From Table | we see that for ConfTT/PhysTT, \Ei,/ii\ > 
il/m)-\ and for CTS, \Eb/^i\ « {l/m)-\ Since gravity 
in general relativity is typically stronger than Newtonian 
gravity, we find again that CTS and ConfTT/PhysTT 
contain too much energy relative to the black hole 
masses, ConfTT/PhysTT having even more than CTS. 

The proper separation between the apparent hori- 
zons £/m is about 4% smaller for CTS than for Con- 
fTT/PhysTT. By Eq. (g) this should lead to a relative 
difference in binding energy of the same order of mag- 
nitude. Since E^/i-i differs by almost a factor of two be- 
tween the different decompositions, the differences in £/m 
play only a minor role. 



^ Although the Misner solution can be obtained analytically, we 
found it more convenient to solve the Hamiltonian constraint 
numerically. The configuration in Table h corresponds to a sep- 



aration /3 = 12 in terms of p^ . 
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TABLE I: Solutions of difFerent decompositions for two black holes at rest. Ham and Mom are the rms residuals of the 
Hamiltonian and momentum constraint, i is the proper separation be twee n the apparent horizons. mConfTT represents the 
modified conformal TT decomposition which is explained in section VC. inv. symm. represents a conformally fiat, time 



symmetric and inversion symmetric solution of the Hamiltonian constraint. 







Ham 






Mom 




Eadm 


Aah 


Mah 


I 


l/m 


EADM/m 


Emprc / Eadm 


Eb/fi 


ConfTT 


9 


X 


10~ 


7 


4 


X 


10" 


7 


2.06486 


57.7369 


1.07175 


8.062 


3.761 


0.9633 


0.2660 


-0.1467 


PhysTT 


9 


X 


10" 


7 


3 


X 


10" 


7 


2.06490 


57.7389 


1.07176 


8.062 


3.761 


0.9633 


0.2660 


-0.1467 


CTS-add 


2 


X 


10" 


6 


4 


X 


10" 


7 


2.08121 


62.3116 


1.11340 


8.039 


3.610 


0.9346 


0.2434 


-0.2615 


CTS-mult 


2 


X 


10" 


6 


5 


X 


10" 


7 


2.05851 


60.8113 


1.09991 


8.080 


3.672 


0.9358 


0.2444 


-0.2569 


mConfTT 


3 


X 


10" 


6 


1 


X 


10" 


6 


2.0827 


62.404 


1.1142 






0.9346 


0.2434 


-0.2617 


inv. symm. 


















4.36387 


284.851 


2.38053 


17.731 


3.724 


0.9166 


0.2285 


-0.3337 



B. Configurations with angular momentum 

Now we consider configurations which are approximat- 
ing two black holes in orbit around each other. The con- 
formal metric is still a superposition of two Kerr Schild 
metrics. The black holes are located along the x-axis 
with a coordinate separation d. For ConfTT/PhysTT, 
we boost the individual holes to some velocity ztwcy along 
the y-axis. For CTS we go to a co-rotating frame with 
an angular frequency fl = flcz- Thus, for each decompo- 
sition we have a two parameter family of solutions, the 
parameters being {d, v) for ConfTT and PhysTT, and 
(d, ft) for CTS. 

By symmetry, these configuration will have an ADM 
angular momentum parallel to the z-axis which we de- 
note by J. In order to compare solutions among each 
other, and against the conformally flat inversion sym- 
metric data sets, we adjust the parameters (c?, v) and 
(d, n), such that each initial-data set has angular mo- 
mentum J/ fim = 2.976 and a proper separation between 
the apparent horizons of Z/m = 4.880. In Ref. ||3^, these 
values were determined to be the angular momentum and 
proper separation of a binary black hole at the innermost 
stable circular orbit. 

Table |l| lists the parameters corresponding to this situ- 
ation as well as results for each initial-data set^. As with 
the configuration with black holes at rest, we find again 
that ConfTT/PhysTT and CTS lead to different ADM- 
energies. Now, EADM/m and Emprc /Eadm differ by 
0.02 and 0.013, respectively, between CTS and Con- 
fTT/PhysTT. However, in contrast to the cases where 
the black holes at rest, now CTS and the inversion sym- 
metric data set have very similar values for EADM/m 
and Emprc / Eadm- 



^ Because of the Lorentz contraction, the apparent horizons for 
ConfTT/PhysTT intersect the sphere with radius 2. In order to 
have the full apparent horizon inside the computational domain, 
the radius of the excised spheres was reduced to 1.9 for these 
data sets. 



C. Reconciling conformal TT and thin sandwich 

We now investigate further the difference between Con- 
fTT/PhysTT and CTS. Since the resulting initial-data 
sets for PhysTT and ConfTT are very similar, we restrict 
our discussion to ConfTT. 



1. Motivation 

The construction of binary black hole data for the Con- 
fTT/PhysTT cases produces an extrinsic curvature that 
almost certainly contains a significant TT component. It 
would be interesting to know how significant this compo- 
nent is to the value of the various physical parameters we 
are comparing. Ideally, we would like to completely elim- 
inate the TT component and see what effect this has on 
the resulting data sets. Unfortunately, this is a difficult, 
if not impossible, task. 

The TT component of a symmetric tensor M*-' is de- 
fined as 

M'rji^ = M'^ ^ {"LYY^ , (85) 

where the vector is obtained by solving an elliptic 
equation of the form 

AlF'^VjM*^'. (86) 

The problem resides in the fact that the meaning of the 
TT component depends of the boundary conditions used 
in solving (|86|). 

For the ConfTT/PhysTT cases we are actually solving 
for a vector that is a linear combination of two compo- 
nents, one that solves an equation of the form of ( |86| ) to 
obtain the TT component of A/*-' and one that solves the 
momentum constraint. But by imposing inner-boundary 
conditions of = 0, we don't specify the boundary con- 
ditions on either part independently. Nor is it clear what 
these boundary conditions should be. Since we cannot 
explicitly construct the TT component of the extrinsic 
curvature, we cannot eliminate it. Although it is not 
ideal, there is an alternative we can consider that does 
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TABLE II: Initial-data sets generated by different decompositions for bi nary black holes with the same angular momentum 
J/yLm and separation i/m. The mConfTT dataset is explained in section VC. It should be compared to CTS-add. 







parameters 


Mah 


Eadm 


J 1 [im 




EADhi/m 


Emprc / Eadm 




ConfTT 


d = 


= 11.899, u = 


= 0.26865 


1.06368 


2.12035 


2.9759 


4.879 


0.9967 


0.2906 


-0.0132 


PhysTT 


d = 


-- 11.899,1; = 


= 0.26865 


1.06369 


2.12037 


2.9757 


4.879 


0.9967 


0.2906 


-0.0132 


CTS-add 


d = 


= 11.860, f7 


= 0.0415 


1.07542 


2.10391 


2.9789 


4.884 


0.9782 


0.2771 


-0.0873 


CTS-mult 


d = 


= 11.750, f7 


= 0.0421 


1.06528 


2.08436 


2.9776 


4.880 


0.9783 


0.2772 


-0.0867 


mConfTT 


d = 


= 11.860, f7 


= 0.0415 


1.0758 


2.1061 


3.011 


4.883 


0.979 


0.278 


-0.085 


inv. symm." 












2.976 


4.880 


0.9774 


0.2766 


-0.09030 



"Data taken from 
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provide some insight into the importance of the initial 
choice of M*^ . 



2. Black holes at rest 



caused by our procedure to determine M'^'' . We force the 
extrinsic curvature of mConfTT into the form Eq. (|8S|). 
This is precisely the form of the extrinsic curvature in 
CTS, Eq. (|2^), even using the same function N and the 
same boundary conditions on the vectors Y"^ and /3*. 



Consider the following numerical experiment for two 
black holes at rest: Given N'P^ from Eq. (^8|), make a 
transverse traceless decomposition by setting 



{lYy 



(87) 



where WjM^j, — and 



iV = iV^ + iV^ - 1. Notice 
that we are decomposing 2NAP^ , not AP^ . Taking the 
divergence of Eq. (p^) one finds 



The decomposition chosen in Eq. (|87| ) is motivated by 
the conformal thin sandwich decomposition. With this 
decomposition we can, in fact, use the shift vector N'^ 
to fix boundary conditions on , just as it was used 



(51b 



0. 



. For 

After 



to fix the boundary conditions in Eqs 
the black holes at rest in this case, we have 
solving Eq. ( p8[ ) for F*, we can construct a new conformal 
extrinsic curvature by 



(89) 



which is similar to what would result if we could elim- 
inate M'jij. from M^-i. Using M"-i in place of M*^ we 
can again solve the conformal TT equations. The result 
of this modified conformal TT decomposition "mCon- 
fTT" is striking: Figure shows that mConfTT gener- 
ates a conformal factor -0 that is very similar to ?/; of CTS. 
mConfTT is also included in Table | where it can be seen 
that the quantities Eadm I'm and Emprc / Eadm differ 
only shghtly between mConfTT and CTS. 

The fact that modification of the extrinsic curvature 
changes the ADM-energy by such a large amount under- 
lines the importance of a careful choice for the extrin- 
sic curvature M'^ in ConfTT/PhysTT. The extremely 
good agreement between CTS and mConfTT is probably 



3. Black holes with angular momentum 

We now apply the modified conformal TT decompo- 
sition to the orbiting configurations of section VB. In 



the corotating frame, the black holes are at rest, and we 
start with 7^- and M'-' of two black holes at rest with co- 
ordinate separation d = 11.860. We now solve Eq. ( |8^ ) 
with 



N ^ Na + Nb -1 

and corotating boundary conditions on 
Eqs. (|5Tb|)-(|5Tdi)]: 



(90) 
[cf 



y = iVJ sphere inside hole A, (91a) 

Y' = N's sphere inside hole B, (91b) 

Y^ = fl X f outer boundary. (91c) 

Na/b and N'^a/b lapse and shift of individual Kerr- 

Schild black holes at rest. M'*-' is again constructed by 
Eq. ( |89| ) and used in solving the conformal TT equations. 

Results from this procedure are included in Table ||. 
Again, mConfTT generates results very close to CTS. 
Eadm I'm changes by 1.8% of the total mass between 
ConfTT and mConfTT, again highlighting the impor- 
tance of the extrinsic curvature. 

D. Dependence on the size of the excised spheres 

The framework presented in this paper requires the 
excision of the singularities at the centers of each hole'^. 



Marronetti and Matzner[Q effectively excised the centers, 
by using "blending functions". 



too, 



15 





FIG. 9: Plots of tp and along the positive a;-axis for Con- 
fTT for different radii Ve^c = 2M, M, 0.5M, 0.2M. The ex- 
cised spheres are centered on the x-axis at x = ±5. The 
position where a line terminates gives r^xc for that line. 



TABLE III: Solutions of solving ConfTT for different radii of 
the excised spheres, Vexc- The results for PhysTT are nearly 
identical. 



Eai 



Aai 







Conformal TT 






2.0 


2.0649 


57.737 


8.062 


0.9633 


3.904 


1.0 


2.0682 


57.825 


8.101 


0.9641 


3.917 


0.5 


2.0808 


58.520 


8.101 


0.9642 


3.893 


0.2 


2.0978 


59.514 


8.093 


0.9640 


3.858 


0.1 


2.1064 


60.025 


8.089 


0.9638 


3.840 



So far we have used r^xc = 2M or r^xc = 1.9Af in or- 
der to impose boundary conditions close to the appar- 
ent horizons, but different choices can be made. Indeed, 
one might expect that the boundary conditions ( ^ ) and 
(|l|) become "better" farther inside the apparent horizon, 
where the metric and extrinsic curvature of that black 
hole dominate the superposed metric 7y and superposed 
extrinsic curvature M'^ . 

In order to test this assumption, we solve the constraint 
equations for two black holes at rest for different radii 
Texc- We find that for all three decompositions, the data 
sets depend strongly on the radius of the excised spheres. 
Figure ^ presents plots of the conformal factor ip and 
for ConfTT with different Texc- There is no clear 
sign of convergence of ^ as rexc 0. For rexc — 0.2A/, 
the conformal factor ip even oscillates close to the ex- 
cised sphere. Table III displays various quantities for 

As 



the ConfTT decomposition for different r^xc- As 
varies between 2. OA/ and O.IA/, the ADM-energy varies 



FIG. 10: Apparent horizons for ConfTT with different radii 
of excised spheres. Results shown are for r^xc. = 2M (long 
dashed line), M (dotted line), 0.5M (short dashed line) and 
0.2M (outer solid line). The inner solid line is a circle with 
radius 2. The insert shows a parametric plot of r{<j)) — 2, 
which emphasizes the differences between the different ap- 
parent horizons. 



between 2.065 and 2.106, whereas the apparent horizon 
area changes by nearly 4%. The apparent horizons move 
around somewhat as r^xc changes. Figure 10 shows the 
location of the apparent horizons for different r^xc- 
For CTS-add (with a = Na + Nb - I), the initial- 



data sets seem to diverge as Vexc is decreased. This has 
to be expected, since this choice for a changes sign if 
the excised spheres become sufficiently small. Changing 
to a = NaNb so that the lapse does not change sign 
reduces this divergent behavior. Von Neumann boundary 
conditions on ip at the excised spheres, 



^ = 0, 
dr ' 



(92) 



lead to an increase in Aah especially for large excised 
spheres. This combination of lapse a and boundary con- 
ditions exhibits the smallest variations in Eadm /'"n; cuts 
through ip, (3^ and through the apparent horizons are 
shown in Figs. |l^ and 12, From the three examined 
combinations of lapse and boundary conditions, the one 
shown behaves best, but there is still no convincing sign 
of convergence. 

Table IV presents ADM-energies and apparent hori- 



zon areas and masses for CTS with different 



and 



different choices of lapse and boundary condition. From 
the unsealed ADM-energy Eadm it is apparent that 
a = Na + Nb — 1 diverges most strongly. Note that 
between all choices of lapse, boundary conditions and 
Texc, the unsealed quantities Eadm, Mah, and £ ex- 
hibit a much broader variation than the scaled quantities 
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X 



FIG. 11: Cuts through i/) and P"" for CTS-mult for different 
radii Vexc- Here a = NaNb and the boundary condition on 
-0 at the excised spheres is dtli/dr = 0. The curves for are 
shifted up by 0.5 for x < 5, and are shifted down by 0.5 for 
a; > 5 to allow for better plotting, dtp/dr approaches zero 
at the inner boundary on scales too small to be seen in this 
figure. 



TABLE IV: Solutions of CTS as a function of radius of ex- 
cised spheres, r^xc- Different choices of the lapse a and bound- 
ary conditions for ip at the excised spheres are explored. 





Eadm 


Aah 


£ 


EADu/m 


£/Eadm 






a = Na + Nb ~ 1, 


■0 = 1 




2.0 


2.0812 


62.312 


8.039 


0.9346 


3.863 


1.0 


2.1846 


68.279 


8.000 


0.9372 


3.662 


0.5 


2.3085 


76.253 


7.925 


0.9371 


3.433 


0.2 


2.5463 


93.534 


7.750 


0.9333 


3.044 


0.1 


2.8543 


118.834 


7.489 


0.9282 


2.624 






a = NANB, V = 


= 1 




2.0 


2.0585 


60.811 


8.080 


0.9358 


3.925 


1.0 


2.1216 


64.080 


8.044 


0.9395 


3.792 


0.5 


2.1696 


66.790 


8.017 


0.9411 


3.695 


0.2 


2.2120 


69.456 


7.991 


0.9409 


3.613 


0.1 


2.2326 


70.809 


7.978 


0.9405 


3.573 






u^NaNb, dip/d 


r = 




2.0 


2.1110 


64.229 


8.085 


0.9337 


3.830 


1.0 


2.1533 


66.128 


8.030 


0.9387 


3.729 


0.5 


2.1794 


67.427 


8.011 


0.9409 


3.676 


0.2 


2.2136 


69.559 


7.990 


0.9409 


3.609 


0.1 


2.2330 


70.836 


7.978 


0.9405 


3.573 



2 
1.5 
1 
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FIG. 12: Apparent horizons for CTS with a = Na Nb and 
inner boundary condition dxp/dr — 0. The different curves 
belong to different reic as explained in Fig. |lo| 



EADAi/m and I/Eadm- 



VI. DISCUSSION 



Our results clearly show that different decompositions 
lead to different initial-data sets, even when seemingly 
similar choices for the freely specifiable pieces are used. 
From Tables | and one sees that Eadm /in changes 
by as much as 0.029 between ConfTT/PhysTT and CTS. 
The difference between ConfTT/PhysTT and the inver- 
sion symmetric data is even larger, 0.047. These numbers 
seem to be small; however, current evolutions of binary 
data usually find the total energy emitted in gravitational 
radiation Eaw/'m to be between 0.01 and 0.03 ||l|, ||, |8|, 
which is the same order of magnitude as the changes in 
Eadm /"tl we find. This means that, in principle, most of 
the energy radiated in these simulations could originate 
from "spurious" energy in the system and not from the 
dynamics of the binary system we are interested in. 

These findings highlight the fact that current binary 
black hole initial data sets are inadequate for the task 
of accurately describing realistic binary systems. We see 
that the choices of the conformal 3-geometry 7^- and the 
freely specifiable portions of the extrinsic curvature, em- 
bedded in M^^ , influence the content of the initial data at 
a significant level. Furthermore, the results suggest that 
small changes in the free data associated with the ex- 
trinsic curvature are more significant than small changes 
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in the choice of 7^ This assertion is supported by the 
fact that Ead Af /™ is consistently larger for the ConfTT 
solutions than for the CTS solutions but the two ap- 
proaches can be made to produce quite consistent results 
by using the modified extrinsic curvature of the mCon- 
fTT method. All of these decompositions use the same 
non-flat conformal metric, but differ in the extrinsic cur- 
vature. On the other hand, results for the conformally- 
flat inversion-symmetric data agree rather well with the 
results from the CTS method when we consider orbiting 
black holes. For black holes at rest, CTS differs from the 
inversion symmetric data, which seems to contradict our 
conclusion. However, this difference is likely due to the 
time-symmetric nature of the inversion symmetric data, 
which is especially adapted to the time-symmetry of the 
particular configuration of "two black holes at rest" . 



Improved binary black hole initial data will require 
choices for the freely specifiable data that are physically 
motivated, rather than chosen for computational conve- 
nience. The same is true for the boundary conditions 
used in solving the constraints. The boundary condi- 
tions used in this paper carry the implicit assumptions 
that the approximate metric and extrinsic curvature are 
correct at the excision boundaries and that the value of 
the single-hole Kerr-Schild shift at the excision bound- 
ary is correct in a multi-hole situation. This is clearly 
not true, but we might hope that the impact of the error 
in this choice would diminish as we decrease the radius of 
the excision boundary. However, our results presented in 
Tables [H and |^ do not support this conjecture. Exam- 
ining the change in Ead m /w as we vary rexc shows only 
a small change, but more importantly, it shows no sign of 
converging as we decrease Vexc- The effects of changing 
rexc are much more significant for £/ Eadm, changing its 
value by as much as 10% in the case of CTS-mult for 
the range of values considered. Furthermore, as with 
the energy, we see no sign of convergence in £/ Eadm as 
rexc decreases. Interestingly, although the solutions show 
no sign of convergence as we shrink the excision radius, 
we do find that the dimensionless quantities Eadm I'm 
and IjEADM do become independent of the choice of the 
inner-boundary condition as Vexc decreases. This can be 



seen in comparing the result in Table IV for the cases us- 
ing ijj = 1 and dip/dr = as inner-boundary conditions. 
Additional tests, not reported in this paper, further sup- 
port this assertion. 



* Follo wing submission of this paper, a preprint by Damour 
et al.|B9|| has appeared that lends support to our idea that 
the extrinsic curvature plays a key role in constructing quasi- 
equilibrium binary black hole initial data. 



VII. CONCLUSION 

Using a new elliptic solver capable of solving the initial- 
value problem of general relativity for any of three differ- 
ent decompositions and any choice for the freely specifi- 
able data, we have examined data sets representing bi- 
nary black hole spacetimes. We find that the choices for 
the freely specifiable data currently in use are inadequate 
for the task of simulating the gravitational radiation pro- 
duced in astrophysically realistic situations. In particu- 
lar, we studied the results of using a superposition of two 
Kerr-Schild black holes to fix the freely specifiable data 
and compared them to the results obtained from confor- 
mally flat initial data. 

Although the new Kerr-Schild based data provide 
a valuable point of comparison, it is not clear that 
the data produced are significantly superior to previous 
conformally-flat data. What is clear is that the choice 
of the freely specifiable data will be very important in 
constructing astrophysically realistic binary black hole 
initial data. Progress will require that these data, and 
the boundary conditions needed to solve the constraints, 
must be chosen based on physical grounds rather than 
computational convenience. 

How can better initial data be achieved and how can 
the quality of initial data be measured? We believe that 
the conformal thin sandwich decomposition will be espe- 
cially useful. Genuine radiative degrees of freedom can- 
not in principle be recognized on a single time slice. The 
conformal thin sandwich method uses in effect two nearby 
surfaces, giving it a potential advantage over other meth- 
ods. Also, it avoids much of the uncertainty related to 
specifying a conformal extrinsic curvature. Moreover, 
the conformal thin sandwich approach is especially well 
suited for the most interesting configurations, a black- 
hole binary in a quasi-equilibrium orbit. In this case 
time derivatives of all quantities are small and the choice 
M*-' = is physically motivated. One should exploit the 
condition of quasi-equilibrium as fully as possible, i.e. 
one should use the conformal thin sandwich approach 
together with the constant K equation, dtK = 0. The 
latter yields another elliptic equation for the lapse which 
removes the arbitrariness inherent in choosing a confor- 
mal lapse a. One will also need more physical boundary 
conditions. Work in this direction was begun in[^ 
and refined in |4^ . 

Ultimately, the gravitational wave content of an initial- 
data set can be determined only by long term evolutions. 
One must compute an initial-data set representing a bi- 
nary black hole in quasi-circular orbit and evolve it. Then 
one must repeat this process with an initial-data set rep- 
resenting the same binary black hole, say, one orbital 
period earlier, and evolve that data set, too. If both 
evolutions lead to the same gravitational waves (modulo 
time offset) then one can be confident that the gravita- 
tional radiation is indeed astrophysically realistic. This 
approach has recently been used for the first time in con- 
junction with conformally flat puncture data, where it 
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proved remarkably successful |3|. 
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